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Abstract. We compare and contrast three different perturbative expansions for the 
quartic anharmonic oscillator wavefunction and apply a modified Borel summation 
technique to determine the energy eigenvalues. In the first two expansions this provides 
the energy eigenvalues directly however in the third method we tune the wavefunctions 
to achieve the correct large x behaviour as first illustrated in [T] . This tuning technique 
allows us to determine the energy eigenvalues up to an arbitrary level of accuracy with 
remarkable efficiency. We give numerical evidence to explain this behaviour. We also 
refine the modified Borel summation technique to improve its accuracy. The main 
sources of error are investigated with reasonable error corrections calculated. 



1. Introduction 

Quantum Field Theory (QFT) is our most successful mathematical framework for 
describing the fundamental laws of physics. Despite this we still have few tools to 
calculate physical quantities in many models. Perturbation theory is probably the 
most important technique available to calculate such quantities however perturbative 
expansions are often divergent within their region of applicability. In addition some 
physical properties are not correctly reflected in perturbative expansions. For example, 
the renormalisation group implies that energy eigenvalues in Yang-Mills theory cannot 
be solved for perturbatively. 

One technique that has been of interest in both the Schrodinger representation of 
quantum field theory and quantum mechanics is a modified type of Borel resummation 
[2] [S] [1] • The aim of this paper is to explore this particular method of resummation. 
It is important to test and develop the accuracy of our technique in a theory with known 
results so that we can confidently apply it to problems in which other techniques fail. 
Since energy eigenvalues are already known for the quantum mechanical anharmonic 
oscillator to a high degree of accuracy (e.g. [1][6] and references therein), we shall apply 
the resummation technique to this problem. That is we look for solutions to 

- — - + px^^ + gx^^ = E-^ (1) 
hm * = (2) 

|3;|— >oo 
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where x is defined along the real axis and we choose units in which h = 2m = 1. We will 
generate three different perturbative expansions of the quartic anharmonic oscillator 
wavefunction as follows: 

• An expansion in the coupling, g (Bender Wu expansion) 

• An expansion in Plank's constant, h (the semi classical expansion) 

• An expansion in powers of x 

In the first two approaches we take a perturbative expansion of the energy eigenvalue 
and directly resum to find an approximate result. In the third method we resum the 
large x behaviour of the wavefunction. This provides a variational technique in which 
the energy is tuned to ensure the correct boundary condition ([2]) is observed. This 
technique is remarkably efficient and provides energy eigenvalues up to an arbitrary level 
of accuracy as first illustrated in pQ. We shall apply the modified Borel resummation 
technique to provide further explanation as to why this method is so efficient. 

The method of resummation we shall employ allows us to extract the small s 
properties from an asymptotic expansion which is only valid for large s. Therefore 
consider an asymptotic expansion of a function f{s) in inverse powers of s: 

/(s) ^ ao + - + ^ + ^ + ^ + ^ + ■ ■ ■ . (3) 

We analytically continue f{s) into the complex s plane, then Cauchy's theorem relates 
the large s to small s = sq behaviour of f{s) via the integral 

L(X)= / ds f{s) (4) 

Jc s — So 

where C is a large circular contour centred on the origin. We only assume f{s) to be 
analytic in the half plane 3?(s — sq) > 0. Any singularity contributions to (jll) in the half 
plane 3?(s — Sq) < are exponentially dampend by A so that lim^^oo L{X) = /(sq). 

We approximate L{X) using a truncated version of the asymptotic expansion 
and expanding the s — sq denominator in powers of so/s truncated to some order P. 
Thus we define 

r p p N 

L,(X)= ds.M'-'-^^J^Y.^ (6) 



c 



p=0 n=0 



n=0 p=0 ^ ^ ' 

where in completing the integral we used the identity ds s"" exp(As) = 27riA"^^/r(n) 
for n < 0. The introduction of the Gamma functions in this series improves the 
convergence of the original asymptotic expansion. We note however that the LfqiX) 
is only a good approximation to Ij{X) within a limited range of A. For sufficiently 
large A the series will be dominated by the highest powers of A and exhibits a rapidly 
increasing or decreasing behaviour depending on the sign of the coefficient. We also 
note the similarity of this method to that of Borel summation. The Borel transform of 
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an asymptotic series results in the introduction of an additional l/n\ factor in each c„. 
The Borel procedure however requires us to analytically continue the Borel transformed 
series before inversion. This is accomplished via techniques such as Pade approximants 
or conformal mapping. The advantage of our technique is that this analytic continuation 
is encoded in the contour integral. 

We take A as large as practicably possible with the constraint that Ljv(A) be a good 
approximation to L{X). This is best achieved by requiring maximal A = Am such that 
LAr(A) differs from Ln^i{X) by no more than a set amount. In this paper we will choose 
Am so that they differ by no more than 10~^ percent. More terms (greater A^) allows 
for a larger A and therefore better dampening of any singularity contributions. 

One method for increasing the singularity dampening for a given is to introduce 
a new parameter, a by replacing f{s) with /(s") and sq with Sg^" in L{X). For a > 1 
the size of the last term in L^r relative to its penultimate term L^-i is reduced due 
to the gamma function in This allows us to take a larger value of A whilst La? (A) 
remains a good approximation to L{X). Increasing a however causes singularities of / 
to be rotated about the origin. For a too large the singularities enter the half plane 
3?(s — So) > at which point they are no longer exponentially suppressed. L{X) will then 
exhibit oscillations resulting from these singularity contributions. We therefore take a 
to be as large as possible but still ensuring Ljv(A) is monotonic as a function of A for 
A < Xm- The technique originally (i.e. without the introduction of a) only worked for 
functions which are analytic in the half plane 3fJ(s — Sq) > 0. With the introduction of a 
we could also consider functions in which /(s) has singularities with 3?(s — sq) > 0. By 
reducing a < 1 we can rotate these singularities back into the half plane 3?(s — Sq) < 
where they become exponentially dampened. 

In all three approaches we will solve the anharmonic oscillator ground state by 
writing = since the ground state of any quantum mechanical system has no 
nodes. In the case of the quartic anharmonic oscillator ([T]), the potential and boundary 
condition are even in x. We therefore expand W in the form W = X^^i'^n^^"- 
Substituting this expansion into the differential equation ([T]) and comparing coefficients 
of the x^" we get relations between the follows 

- 2ai = E, - 12a2 - 4a^ + p = 0, - SOag - 16aia2 + g = (7) 

and for n > 3 

a„+i = - I ^4m{n - m + l)a„a„_m+i | /{2{n + l){2n + 1)). (8) 

\m=l / 

The equation ([H]) allows us to find a„+i in terms of the with m < n. We could 
therefore solve all of the a„ in terms of oi, 02 and 03. In turn these first three coefficients 
are determined by the physical parameters E, p and g via ([7]). However for a given p 
and g only specific values of E allow the solution to satisfy the boundary condition ([2]). 
In sections |2l [3] and H] we propose different methods for eliminating this final degree of 
freedom. 
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2. Bender Wu Expansion 

Bender and Wu [7] showed how to construct the ground state energy by summing all 
connected Feynman diagrams with no external legs. In general a Feynman diagram with 
2n external legs has at least n — 1 vertices. So we ensure that an x^" term is at least 
order n — 1 in the coupling by making the expansion a„ = Ylm=n-i ^n,mg"^- This is a 
similar approach to another method outlined by Bender and Wu in the same paper. We 
substitute this coupling expansion into the above relations between the a„ and compare 
coefficients of g"'. The first coefficient is given by a^Q = 1/4 which requires a choice of 
sign for aifi. In keeping with the Bender Wu methods |7]| we choose a negative sign. In 
[1] we showed that this sign choice is required to ensure the correct boundary condition 

Having made the choice for ai^o, the remaining coefficients are uniquely determined. 
We ffist find each an,n-i by looking at the g'^~^ coefficient in ([7]) and dH]). The g"' 
coefficients give an,n then the f?""*"^ coefficients give the a„,„+i etc. At each stage we are 
substituting in the previous solutions. Eventually we can find each a„ up to any order 
in g. Of course the coefficients will also depend on p however for the purpose of this 
section we set p = 1 without loss of generality since the eigenvalues for arbitrary p may 
be reproduced via a form of Symanzik scaling [S] . 

We now have a solution for the ground state wavefunction of ([T]) as the exponential 
of a power series in x and g. The energy eigenvalue is computed from ai using ([7]) as a 
power series in g. This gives the well known [7] [9] Rayleigh Schrodinger perturbation 
expansion for E{g), 

^ , 3 21 2 333 o 30885 . 916731 , 65518401 ^ 

E =l + -g g^ H g^ g^ H g^ g^ + . . . .(9) 

4^ 16 64 ^ 1024 ^ 4096 32768 ^ ^ 

This expansion has a zero radius of convergence as can be seen from the asymptotic 
form of the ai,„ coefficients at large n as given by Bender Wu [7]. It has been used to 
generate some energy eigenvalues via a Borel resummed Fade approximants technique 
[TT][T2] although many other techniques have been used to find the energy eigenvalues 
more accurately and efficiently e.g. [T] p] [H] [15] pT] . 

We will apply our resummation method to E{g) in an attempt to get meaningful 
results for non zero values of g. To do this we analytically continue g in the complex 
s = 1/g plane. We write Sq = 1/g and c„ = — 2ai_„ then apply the contour integral 
technique so that Ln{\) in ([5]) approximates the energy. 

The results generated via our resummation method are listed in figure [1] and 
compared to the results as generated by the method [T]. We label the eigenvalues 
generated by [I] E^^^^ since they are accurate to within the number of significant figures 
expressed. The seemingly strange choices for g that we use becomes more natural in 
the semi classical expansion as outlined in the next section. We use the same values of 
g in both sections to allow comparison. 

The results for small g are quite impressive with errors in the region of 10^^ to 10~^ 
percent however for larger g the results are less impressive with the error approximately 
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Figure 1. Results for rcsummatfon of E in tlie coupling. 



23.8% for the largest value of g. 

Figure [2^ (a small g example) shows the expected behaviour of L2g{X) and L3o{X) 
with P = 50. For sufficiently large A the two curves become a good approximation 
to L{X) and we see a flattening of the curve. For A sufficiently large we notice an 
appreciable divergence of the two curves. Figure [2)d where g is relatively large has 
somewhat different behaviour. Here we notice that the curves have not started to 
flatten before they appreciably diverge. For these larger g we need more terms (greater 
A^) in the expansion so that we may consider larger A where the curve starts to flatten. 
Additionally we note that this curve exhibits oscillatory behaviour although remains 
monotonic. When we introduced the requirement for the curve to be monotonic we 
assumed that a singularity contribution would consist of an exponentially weighted 
sinusoidal correction to a flat curve. In this case we do not have sufficient terms to 
consider A in the region where it becomes flat but instead are considering a region of 
the curve where it is still appreciably increasing. If we considered more terms with this 
value of a we may well flnd our monotonic condition is violated. Given the behaviour 
observed in flgure [2b it is not surprising that we have such a large error. 

In the next subsection we show thal|j E ~ g^^^. Therefore as (? ^ 00 or equivalently 
So ^ we notice that E has a singularity. It is this singularity that is causing difficult 
in the resummation process since we need larger A for large g to ensure it is damped 
sufficiently. This in turn requires a larger A^. The problems no doubt could be solved 
if we took a sufficient number of terms in the expansion of E however this would be at 
the expense of greater computing resources. In the next subsection we outline a more 

I This can be seen by Symanzik scaling of the Hamiltonian [8] . The infinite coupling limit has been 
calculated in \T^ . 
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a) g = 0.05783 



b) ^ = 206.09853 



Figure 2. Plot of i3o(A) and £29 (A) with a = au 

efficient method in which the singularity contribution is absent. We will therefore be 
able to calculate the infinite coupling limit. This is something that cannot be done via 
resummation of the Bender Wu expansion. 



3. Semi Classical Expansion 

In this section we resum a semi classical expansion in h for the a„ coefficients in W . We 
write = e^l^ and consider the modified differential equation 



+ (60 + feix' + h^x") ^ = 0. 



(10) 



The original problem ([T]) in which = m = 1 is then recovered via a rescaling x ^ cx 
provided 



(? 



h^' —-^2' '-"'h^- ^^^^ 

We substitute \1/ into the h dependent differential equation to generate relations between 
the On- The equivalent expressions to (JTj) are 

bi = 4al + 12a2h, 



bo = 2aih, 
and the new version of (IHl) is 



62 = 16aia2 + SOa^h 



2h{n + l){2n + l)a„+i + ^ 4m(n - m + l)am,a„_m+i = 0. 



(12) 



(13) 



m=l 



The advantage of performing this rescaling is that we now have some freedom to choose 
ai and 02. We will restrict the choice however by requiring ^ > with both c and h 
real. This will allow us to choose Oi and 02 up to a sign. We will choose Oi = —1/2 and 
02 = —1/8. We showed that this was the appropriate sign choice in [1] to ensure that 
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Figure 3. Ratio of coefficients of h" in 62 



the boundary condition 
to 

1 

E = 



is satisfied. With these choices we can reduce ffTTll and ffT2 



(14) 



We now assume an expansion a„ = ^ 



(7 fl" 



for each n > 3. This is substituted 



into ( 1131) and coefficients of h"^ compared for each n. We first compare coefficients of 
h'^ to get the Onfl then the coefficients of give the a^^i etc. At each stage previous 
results are substituted into the new equation. By continuing this process sufficiently 
many times we can find each a„ to any order required. A simple program can therefore 
be created to calculate 62 (by using the H expansion of 03) as an expansion in h. The 
first few orders are 



1- 



8 32 



2555 



69545 



4849705 . 



202337485 



h'^+ h^+- ■ ■ (15) 

512 2048 16384 65536 ^ ^ 

The expansion of 62 in h is an alternating sign series so we plot the ratio of successive 
coefficients of in 62 in such a way that we are dividing by the preceding term and 
removing the minus sign. We illustrate this in figure [3] for the first 100 coefficients and 
note the approximate linear behaviour for large orders. This suggests the asymptotic 
behaviour of the coefficients in 62 have the form (— l)"+^A;"r(n + /C2) where ki and /c2 
are real constants. Such asymptotic expansions have a zero radius of convergence as 
demonstrated via the alternating sign series test. 

We shall therefore resum (fT5l) for a particular sq = 1/h from which E and g can 
be calculated. The in L]\r{X) now correspond to the coefficients of h"' in the 62 
expansion. We note that ^ = corresponds to g = and E = 1 whilst /i — > 2/3 
corresponds to (7 — 00. We actually find that h G [0,2/3) corresponds to g & [0, 00) 
which can be confirmed by the results of pj . We therefore calculate some couplings and 
energy eigenvalues within this range. The results are given in figure H] again compared 
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Figure 4. Results of resummation in the semi classical expansion 



to results, gbest (accurate to the stated number of significant figures) generated from pLj. 
The energy eigenvalues produced are exact however it is the coupling that we are trying 
to approximate via the resummation process. We note that errors in the coupling are 
in the order of 10""*%. Most importantly however the error remains within this order 
of magnitude for the full spectrum of h in contrast to resummation in the coupling. We 
can attribute this success to the fact that h G [0,2/3) as opposed to g E [0, oo) and 
the resummation process being most effective for small h or g. Also we encountered 
difficulties in resumming the coupling expansion for large g due to the singularity at the 
origin in the s plane. This problem has been removed in the semi classical expansion. 

For higher h we expect the error to be greater since the contribution from 
singularities increases. It is interesting to note however that whilst this is true for 
the larger values of h, the highest error is when h = 0.05. We attribute this to an 
insufficiently large value of P. We plot L^oW and L29{\) in figure with P = 50 and 
note the decaying behaviour of the curves. With P = 100 say we recover the expected 
behaviour of a flattening curve followed by the divergence of the two curves. One curve 
increases whilst one decreases from the point of divergence as illustrated in flgure [5b. 
This is because the expansion of (s — s^")"^ in powers of s\l°' /s is only valid for large 
s. Whilst this is a valid assumption given the contour of integration, the series does 
require more terms to achieve a suitable level of approximation when sq becomes larger 
or equivalently h becomes smaller. The error as a result of truncation in this expansion 
is systematic hence the decaying nature of the curve for larger values of A. This effect 
becomes more pronounced for larger a. It is possible that a becomes sufficiently large 
to cause this decaying behaviour before singularity contributions becomes significant. 
However in this case the curve will still fail the monotonic condition. We resultingly 
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Figure 5. Plot of L3o(A) and L29(A) with h = 0.05 



take both a smaller a and a smaller A and therefore get greater singularity contributions 
than if we had taken a larger P. Despite the improvement expected with larger P we 
are still able to extract good approximations with P = 50. 

The case ^ — > 2/3 is particularly interesting because this corresponds to the infinite 
coupling limit. With E^o = (3/62(2/3)/2)^/^ this limit corresponds to E ^ E^g^/^. 
Using our resummation method we find E^ = 1.0603632150 compared to the value 
Eoo = 1.06036209 which is exact to the stated number of significant figures. That is an 
error of approximately 1.06 x 10~^%. Parisi [19] was also able to calculate this limit as 
E^ = 1.06038. 

The method outlined so far for the ground state can be generalised to find the energy 
and wavefunctions of the excited states. We write the qth excited state \l/ = P^^'o with 
energy E = Eg + Eq. Now consider 

- /i'^ + (bo + h + bix^ + b2X^)^ = (16) 
which using (fTUl) reduces to 

+ (^^) 

We scale a; ca; to recover the AHO oscillator ([T]) provided that in addition to fll4p we 
also have 

b^c^ bz b^ 
Eg = — = -jEo = . (18) 

We now need to find 63 and Ug using (fTTI) which is easily done using a similar 
approach to that described for the ground state. That is expand 

00 00 

U= ^^c„,^n'"x", (19) 

71=0 m=0 
00 

b3= J2b3,mh'^, (20) 

771=0 
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and substitute into (|T71) . By comparing coefficients of h'^ and x"' we get 

k 

- X] b3,iCn,k-i + {n + 2){n + l)c„+2,fc-2 

. ^, ■ (21) 

{n+l)/2 ^ ' 

+ 4 ^ ^ m{n + 2 - 2m)ara,iCn+2-2ra,k-i-l = 
m=l i=0 

for k > 1. The /c = case reveals 63 = or co,n = for all n. Clearly the former is 
required. The k = 1 case is given by 

(n+l)/2 

- h,iCn,o + 4/1 ^ m{n + 2- 2m)am,oCn+2-2m,o = 0. (22) 

We consider this equation for each n starting with n = and find at each stage that 
Cn,o = until we reach a point where 

63^1 = 4nai^o = — 2n. (23) 

Clearly 63^1 = is required for the ground state, g = 0. Note also that a Taylor expansion 
of (ITSjl would reveal 63,1 = -Eg,o- This corresponds to the gth harmonic oscillator with 
energy eigenvalues E = 2q + l or Eg^ = 2q. So we have Cnfl = for n < g and 63^1 = 2q. 

We can fix one of the coefficients in Uq by a choice of normalisation and therefore 
set Cg^o = 1 with the remaining Cq^n = 0. Now consider the k = 2 case for each n then 
k = 3 etc. We see that at each stage (J2T]) can be solved for Cn,k-i provided n 7^ g or if 
n = qwe get 63,^. 

By solving a series of linear equations we have managed to determine 63 and Uq 
and hence found the excited wavefunctions and energies up to any order required. We 
can then apply our resummation method to evaluate the series expansion. 

One advantage of this method is that we can calculate the energy eigenvalues for 
an infinite coupling. We have already seen that this is given when h = 2/3. By using 
( ITSj) and ( fT4|) we can write 

So in the infinite coupling limit, h 2/3 we have Eq = Eq ^oQ^^^ where 

^,,00 = -^63(2/3)^0,00. (25) 

Applying our resummation method we can calculate 63 and hence Eq^^o- A table listing 
the ffist ten 63 is given in figure 



4. Tuning the Boundary Condition 

In this section we use our resummation technique to try and explain why the method 
for generating energy eigenvalues in p] is so efficient. In this approach we examine the 
large x behaviour o{W = 'Y^^=i CLn^^ by applying the resummation technique already 



A Modified Borel Summation Technique 



11 



n 






^3 


Y^best 


Error (%\ 


1 


3 3554 


1 7 9902 


-1 722251 

X> 1 ^^^wX 


-1 722249 


0001 

w .www X 


2 


3 2057 


1 5 4576 


-4 020691 


-4 020850 


0040 

w . w w^w 


3 


3.1254 


13 8688 


-6 654020 


-6 654572 


0083 

w • w WtJ o 


4 


3 0734 


1 2 7075 


-9 555955 


-9 557404 


0152 


5 


3.0341 


11.7759 


-12 681877 


-12 686239 


0.0344 




3 0046 


1 1 0209 


-16 000874 


-16 012209 


0708 

W . W 1 


7 


2.9831 


10.4064 


-19 489826 


-19 514237 


0.1251 


8 


2.9675 


9.8993 


-23.130233 


-23.176133 


0.1980 


9 


2.9560 


9.4721 


-26.906049 


-26.984993 


0.2926 


10 


2.9473 


9.1068 


-30.803510 


-30.930247 


0.4098 



Figure 6. Infinite coupling values of 63 



outlined. The differential equation implies large x asymptotic behaviour of the form 
W ~ ±^/gx^/3 and ([2]) then requires us to take the solutions that have a negative sign. 
We again employ Cauchy's theorem defining 

so that when the boundary condition is satisfied we have limA^oo^l-^) = ^\^/^^ 
provided W has no singularities in the right half s = 1/x plane. 

In the case of the quartic anharmonic oscillator, W exhibits analytic behaviour in 
the whole complex s plane however with higher order polynomial potentials this is not 
the case. A prescription in which we take W{s) and 1/A^ 1/A'^" in ( 126|) was 

introduced. By reducing a < 1 singularities are rotated into the half plane 3fJ(s) < 
in which case they become exponentially suppressed. Since we will restrict ourselves to 
the quartic oscillator in this paper we shall take a = 1. 

A rescaling x ^ cx {c ^ W) in the differential equation allows us to fix the ratio 
k = ai/a2 = ±4. We make a choice for 02 and sign of k then solve the remaining 
in terms of 03 via the recurrence relation ([H]). We vary 03 until we observe the correct 
boundary condition in T(A). As with the previous direct resummation methods we use 
a truncated expansion of W to approximate T(A) with T/v(A). 

When T/v(A) is plotted for 03 too small to satisfy the boundary condition, we find a 
rapidly decreasing curve. If 03 is too large then we find a rapidly increasing curve. We 
are therefore able to tune 03 by adjusting upper and lower bounds to find an interval 
within which 03 lies. We start with a modest value of N and tune 03 until T/v(A) 
produces a curve flattening for the larger values of A. The range of A considered should 
be chosen to ensure T/v is a good approximation to T. That is up to a value of A at 
which T/v and T/v-i appreciably diverge. To further increase the level of accuracy we 
need to consider whether the curve is rapidly increasing or decreasing for larger values 
of A. We therefore need to increase N. So N effectively determines the level of accuracy. 
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With the sign choice of k and value of 02, 03 is determined up to an arbitrary level 
of accuracy. The physical parameters are then found via 

„ —2ka2 1 , 

E = -gs 27 

(16fcai + 3003)3 

p=^^f^i±^gl. (28) 
(16fcai + 30a3)i 

The procedure can be generalised to find the excited energy eigenvalues. As with 
the semi classical expansion we write the gth excited state \E'g = Pq'^o and energy 
E = Eo + Eg. Then Pg in expanded in powers of x, Pg = Yl'^=o CnX"'. The even powered 
coefficients are set to zero for q odd and the odd powered coefficients to zero for q even. 
Again we used a scaling of the differential equation x ^ cx and substitute in "^g. We 
compare coefficients of x" to get a recurrence relation for the c^. This also determines 
Eg in terms of the c„. We set the first non zero coefficient in the expansion (cq or ci) to 
unity and use the recurrence relation to solve the c„ in terms of C2 or C3. We now define 
U via 



1 r . e 



As 



and use the truncated expansion for Pg to approximate U{\) with Un{\). 

We vary either C2 or C3 as before and again observe the rapidly increasing or 
decreasing behaviour for solutions which do not correspond to eigenstates. For example 
in the odd case we vary r = — C3. There will be multiple values of r which correspond 
to different levels of excitation. We label these r„ such that r„_|_i > r„. With r < Ti. 
U{X) corresponds to a rapidly increasing curve and with ri < r < r2 we get a rapidly 
decreasing curve. With T2 < t < T3 U (A) switches back to a rapidly increasing curve. 
This switch between rapidly increasing and rapidly decreasing occurs whenever we 
increase r past one of the r^. The level of truncation determines the level of accuracy 
as before. 

Although the ground state was analytic in the whole complex s plane for the quartic 
oscillator we found that the excited states were not. Therefore when we tune the curve, 
we get oscillations instead of a flattening behaviour. We could in principal use our 
technique of decreasing a < 1 to move the poles into the left half plane where they are 
suppressed by the exponential factor. This reproduces the flattening behaviour of the 
T/v(A) curve. We find however that it suffices to tune for the oscillations. 

One of the advantages of this approach compared to direct resummation of a 
perturbative expansion in hor g is that we are able to produce solutions in which p < 0. 
These solutions are dominated by instanton effects which are non perturbative in both 
h and g. Also this process allows us to determine an energy eigenvalue to an arbitrary 
level of accuracy. It is the switch between the rapidly increasing or decreasing behaviour 
of Tjv or f/jv that makes this method work so well. In P we postulated a reason for 
this and in this section we use our resummation techniques to further investigate this 
behaviour. 
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4.1. The Ground State - Zeros, Poles and Cuts 

If we solve the differential equation ([T]) for its large positive x behaviour without the 
boundary condition ([2]) we find two asymptotic solutions, W = ±-^x^/3. The general 
large x solution to ([T]) is therefore of the form 

= exp (-^xA +Aexp ( ^xA • (30) 



The boundary condition ([21) however requires us to take A = 0. in general has zeros 
at specific points in the complex x plane depending on A. The only exception to this 
is when A = at which point we only asymptotically approach zero as x 00. For 
A < these zeros lie along the real axis whereas for A > they lie off the real axis, 
somewhere in the complex x plane. There will also be contours in the complex x plane 
along which is purely real and negative. The zeros and negative regions in will be 
manifested as cuts and poles in log\l/;. 

It is clear that any pole or cut in the right half complex x plane would result in 
an exponentially increasing or decreasing Lj^{X). Oscillations will occur due to a pole 
lying off the real axis however we took modest to ensure we only see the beginning 
of the oscillation and hence the appearance of a rapidly increasing or decreasing curve. 
We hypothesised in pQ that this is what is being observed in L]\r{\). In this section we 
present numerical evidence to support this. We are essentially tuning the solution until 
A = 0. For A ^ the resummation is conveniently spoilt in such a way that we are able 
to refine the solution. In this section we will employ Cauchy's theorem to determine 
the location of zeros in the solutions to the differential equation as constructed in the 
previous section. We will then be able to observe the dependence of these zeros on 03. 
Whilst we have so far only presented an argument based on the large x behaviour of \l/ 
we will find that the location of zeros in do indeed correspond approximately to the 
location of zeros in the full solution. 

We note that since the coefficients a„ in W (and a Taylor expanded Wi) are all real, 
any cut or pole in the upper half s plane should necessarily be mirrored in the lower 
s plane. Let us momentarily assume that in the right half s plane W{s) only contains 
poles located at Sp and s*. We then construct 

^(A) = / ds- W(s) 

^ ' 2m Jc s-so ^ ' (31) 

~ W{so) + ccos(As^ + z/)e^('?-'«) 

where c, z/ are real numbers and Sp is split into real and imaginary parts, s^ + is^. If we 
had included a cut contribution instead of a pole then we would expect contributions 
from the whole cut with each portion of the cut weighted by the exponential factor. In 
the zero lies on the right end of a region of negativity corresponding to a potential 
cut in Wi and therefore is the dominant singularity contribution in our contour integral 
ifW = Wi. So whilst W{s) does not necessarily have just one pole in the positive right 
quadrant we will model the cuts by a pole representing a weighted average. Since the 
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Figure 7. Location of zeros 

zero is the dominant contributor we will find that S{X) at least in the case of Wi allows 
us to reasonably approximate the location of our zero. We will apply this model to our 
original ^ and predict zeros close to the zeros of ^>i. As usual we use the truncated 
expansion of W to approximate S{X) by Sn{X)- 

With the above construction we are in a position to numerically determine the 
location of the conjugate poles. Initially we will assume that Sp 7^ and plot Sn{X) for 
a given sq. We expect to see oscillations in the curve either growing (s^ > sq), decaying 
{Sp < Sq) or fixed in amplitude. We tune Sq until we achieve the latter at which point 
we have = sq. With this sq we can determine the frequency of the oscillations in 
Sn{X), this gives us s^. 

If we do not observe the oscillatory behaviour in Sn{X) even with A^, P sufficiently 
large then it may be that Sp = 0. If this is the case then we use Sn{X) to approximate 
W{so) for 60 > Sp. As So becomes closer to Sp however the effect of our pole becomes 
more significant since the exponential dampening factor becomes reduced. We must 
therefore take larger A which in return requires larger N. Alternatively we can try 
increasing a > 1. This causes the final term in Sn to be reduced in comparison to the 
penultimate term due to the F function. This again allows a larger A and hence greater 
exponential dampening. Ultimately though we are unable to extract W{so) from Sn{X) 
for So ^ ^-iid indeed the process becomes more difficult as So approaches Sp. We 
can however determine Sp by approximating ^ to a linear behaviour in this region. We 
justify this by solving the differential equation in a region where W and W" are more 
dominant than any potential term. Numerically this linear behaviour appears to work 
well at least within the region that we were realistically able to plot. 

For 03 too small we find that the zeros lie along the real s axis with zeros approaching 
s — (i.e. X — 00) and indeed attaining this value as 03 becomes appropriate for the 
correct boundary condition to be satisfied. For 03 too large we find that the location 
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of zeros lie off the real axis and into the complex plane. They appear to lie on some 
countour which approaches s = as 03 approaches its correct value. We plot the location 
of some of these complex zeros of \1/ in figure [^. We have restricted this plot to one 
quadrent only however the plot would necessarily be mirrored into all four quadrants. 
The solid line represents the location of zeros of as constructed at the beginning of 
this section. 

We note that the two sets of zeros are not in exact agreement although there does 
appear to be a correlation. We attribute the differences largely due to being a small s 
approximation and we therefore expect the approximation to improve as |s| ^0 which 
appears to be happening in the plot. Unfortunately it becomes increasing difficult 
to numerically calculate zeros of \E' as we get closer to the origin since the period of 
oscillation becomes very large. In order to measure this period we need to be able to 
plot at least the first few oscillations. This requires an increasingly large A which in 
turn increasingly requires a large number of terms. 

We also question how accurately our numerical technique can accurately determine 
the location of zeros. As an example we calculate analytically the location of the zeros 
in \1/ = exp(x^) + exp(— x^). Applying the resummation process to log(\l/t) expanded in 
positive powers of x we are able to numerically calculate the location of the zero. We 
find analytically that a zero exists at 0.745 + 0.43i whereas numerically we find it exists 
at 0.738 + 0.429z. Given the scale of the plot this error is relatively insignificant. 



4-. 2. Excited States - Large x behaviour 

We substitute = Pq'^o into the differential equation ([T]) to produce a new differential 
equation satisfied by Pq, 

d^P ^^"^ + E P - (32) 

dx"^ dx dx 

Taking W ~ —^/gx^/S, this differential equation exhibits two types of large x asymptotic 
solution 

We necessarily choose the first of these to ensure the boundary condition ([2]) is satisfied. 
When the boundary condition is satisfied produces a fiat curve (having adjusted a 

appropriately) which would be expected if we have correctly chosen the large asymptotic 
behaviour. If the second type of asymptotic behaviour is chosen then we would expect 
to observe a correction to T/v(A) due to a singularity at s = 0. We hypothesised in [T] 
that when C2 or C3 does not correspond to an energy eigenstate then we are observing 
the second type of asymptotic large x behaviour. 

Using the resummation technique already outlined, the prefactors corresponding to 
different levels of excitation may be plotted. The prefactors corresponding to the 1st, 

§ We have actually scaled the complex plane so that y/16aj7i2~\-'SQa^s'^ / 3 for simplicity of 

calculation in figure [7] 
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Figure 8. Prefactors corresponding to tlie first tliree excitations 

2nd and 3rd excited states are shown in figure [H] both on a large and small x scale. Since 
these functions are either odd or even we restrict the domain to a; > 0. As would be 
expected the prefactor corresponding to the first excited state has one zero located at 
the origin. The 2nd excited state has two zeros and the 3rd excited state has 3 zeros 
when considered along the whole of the real axis. The tuning procedure did not provide 
a method for determining which level of excitation we have, just that we had determined 
an energy eigenstate. Using out resummation method to plot the prefactor we are able 
to confirm which energy eigenstate has been found by counting the number of nodes. 
The technique in section [H] provides an alternative method. We could use resummation 
in a h expansion to approximate C2 or C3 and then use the tuning method to determine 
the value more accurately. 

We note that for large x these prefactors asymptotically approach some constant 
value as predicted. For example the prefactor corresponding to some energies either side 
of the first excited state energy level are shown in figure For ti < r < T2 we see that 
P quickly becomes very large. This is similarly true for r < ri however this time P 
becomes large in the negative direction. This appears not to be the result of a pole or cut 
otherwise for large x the resummation technique would exhibit singularity contributions. 
Instead we find the prefactor is simply resumming to large values and these values are 
increasing in size rapidly. We test the hypothesis that we are observing prefactors with 
large x behaviour determined by the second type of asymptotic behaviour in (|33i) by 
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Figure 9. P for values of ci either side of d = -0.1458432840772 
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Figure 10. Q for values of a eitlrer side of n w -0.1458432840772 



plotting 



Q = log\P\--^x\ 



(34) 



If the hypothesis is correct then we expect Q to approach a constant for large x. We 
plot these graphs for those same values of Ci as before and display them in figure [TUl 
As predicted these plots do flatten out for large x at least on the scale of exp(x'^). This 
numerically supports our hypothesis. 
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5. A Shifted Expansion 

Let us momentarily remove a from L{X) by setting a = 1. We were unable to directly 
evaluate L{X) even having replaced f{s) with the truncated asymptotic expansion. This 
was due to the s — sq denominator in the integrand which we expanded to give a sum 
of integrals of the form 

ds — = T—— (35) 
c r(n) ^ ^ 

which are more easily evaluated. Unfortunately the expansion of this denominator 
required summing over two variables in L]\f{X). We truncated the expansion of the 
(s — so)~^ term to a relatively high order to avoid complications arising from a truncated 
form of this expansion. This effectively meant we had to compute a series involving a 
sum of NP terms. With our chosen and P this amounted to 1500 terms. 

Instead we shall briefly investigate the possibility of shifting s ^ s + sq in L{\) and 
then using the truncated expansion of f{s) to approximate L{X). That is we write 

oo oo _ 

/(, +so) = y = y^ (36) 

(s + So " ^ s« ' ^ 

71=0 n=0 

by expanding the (s + so)~" terms in powers of sq/s. Our new coeflicients Cn would 
then be dependent on sq so a separate expansion would be required for each sq we try 
to evaluate. 

We then truncate the new expansion and substitute into -^^(A). Having done this 
we reinsert the parameter a with the substitution s — > s". 

^ \ctn 
n=0 ^ ^ 

The advantage of this type of expansion is that we no longer need to worry about a 
double summation truncated to order and P. This reduces the computational time 
to evaluate the series but also prevents the problems encountered in [3] as a result of 
truncation in P. 

We should note however that the a parameter in the original Lat ([5]) is different 
to the a in flH7|) . In the original formulation a caused singularities to be rotated about 
the origin. This is still the case in the new formulation however we now have a different 
origin since s has been shifted. It will depend on the location of singularities as to which 
method will allow a larger a and therefore better dampening of pole contributions. 

We complete our resummation process using the shifted h expansion and present 
the results in flgure [TTJ We note that the errors are considerably larger than in the 
previous method. So whilst this method does have some advantages it is clearly less 
eflicient for the semi classical expansion. We attribute this to the different geometry 
involved when rotating poles by varying a. It is clear that in the shifted method we 
have smaller values of aM and therefore Am which results in less dampening of any 
pole contributions. Although this method is inferior for the expansion of 62 we should 
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Figure 11. Results of rcsumniation in the shifted semi classical expansion 



note that depending on the location of singularities it may prove to be better in other 
expansions and therefore should not be ignored. 

We could question however whether it is fair to compare this method with = 30 
to the previous method which effectively had NP = 1500 terms. In the shifted example 
we could afford to take more terms given the reduced computational power required 
to perform the resummation. We note however that in some perturbative expansions 
calculating an expansion to higher orders can often be the limiting factor. 

6. Error Estimates 

There are two main sources of error in our prescription for evaluating /(sq). The first 
is due to Ln{X) only being an approximation to L{\). We require Ln{X) to differ from 
L7v-i(A) by no more than say a percent. Therefore we can think of a as being an error. 
Smaller a requires evaluating L{X) for a smaller A which has the unwanted side effect of 
reducing the exponential dampending of any singularities of f{s) in the left half plane. 
This is the second source of error. We shall outline a method for estimating the error 
in L{X) as a result of singularity contributions. 

Consider a pole of /(s) located at Sp. We will only be working with asymptotic 
expansions for which all coefficients, a„ are real. This implies that for a pole to exist 
at Sp we must necessarily have a pole at the conjugate location s*. We shall for now 
assume that f{s) has just one pair of poles in the left half s plane but is analytic on the 
remainder of the complex plane. The pair of conjugate poles gives a contribution to (jll) 
of the form 

ccos(A2/p + z/)e^(^^-'o^") (38) 
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where s^" is split into real and imaginary parts s^" = Xp+iyp and c, u are real numbers. 
In the large A limit this contribution approaches zero provided a is not taken so large 
as to rotate the poles into the right half plane. For a general a the correction fl38l) looks 
like an oscillating curve either growing or decreasing in amplitude. For Xp > s^" the 
oscillations are growing and for Xp < s^^" the oscillations are being damped. We will 
use this to tune a until the oscillations are fixed in amplitude, say a = ay at which 
point Xp = sl^". Having determined Xp we determine i/p by calculating the period of 
oscillations. This is best achieved by considering the A derivative of Lf^{X) and looking 
for zeros with a = ar- The phase u and amplitude c are then easily calculated by fixing 
the location of zeros in the derivative of Ljv(A) and also ensuring the correct amplitude 
of the sinusoidal curve. Now with a = the pole contribution is 

Pt + CO = c cos(Xyp + u) (39) 

where pt is the residue of /(s°^) at s^""^ and cc denotes the complex conjugate. So we 
have determined Pt numerically in terms of the parameters c, tjp and u. 

We want to evaluate L{X) (or Lj\f{X)) when a = aM, the maximal value of a for 
which the curve is still monotonic. With this value of a the pole contribution is 



A( Sc~S, 







Pm "" ' 1/... + (40) 



1 / \ a-r/o 



+ ^Vp (41) 



where 



and pm is the residue of /(s"*^) at Sc- We can relate px and pm via 

and therefore can calculate the pole contribution from fl4Up . We think of a as being the 
error in approximating L{X) with Ln{X) which in this case was chosen to be a = 10~^% 
and fHUl) as being the correction to L{X) as an approcimation to /(sq) as a result of a 
pole pair contribution. 

In general f{s) will not have just one pair of poles but a number of poles and 
cuts. The contributions from these singularities will be exponentially weighted with 
singularities lying furthest to the right being most dominant. Due to this exponential 
weighting a cut contribution will look like a pole dominated by the right most end of 
the cut. We can therefore apply our procedure assuming the existence of just one pole 
pair and then interprate the correction (1401) as being the dominant or leading order 
correction. We note however that the right most end of a cut can actually change as 
it is rotated by increasing a. The left and right ends of a cut may actually switch if a 
becomes too large. In this case the error estimate may not be the most dominant but 
will be a more minor correction. 

In order to apply our technique we will require sufficient terms to ensure at least 
two peaks or troughs for a = ut- These must exist within a region where Lm{X) is a 
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Figure 12. Corrected couplings due to dominant singularity contributions 



good approximation to L[X). Clearly = 30 in the Bender Wu expansion is insufficient 
for large couplings as Ln{\) did not approximate L{X) sufficiently well. We also note 
that with the semi classical expansion, the errors where of the order 10~^%. So in this 
case the dominant error is provided by a. In fact the exact values lie between L30 and 
L29 in this expansion so a is actually an error bound. It would however be pointless to 
calculate the error correction due to singularities without ffist reducing a (or A^). We 
will instead apply this technique to the shifted semi classical expansion (section [5]). In 
this technique we found the actual error to be much greater than cr. 

We subtract the dominant singularity correction and list the results {gcorr) in figure 
[T2I The percentage errors after resummation and percentage errors having subtracted 
the dominant singularity contribution are listed. The relative error is the ratio of these 
two errors and indicates that the percentage error has improved by a factor of between 
90 and 312. We expect the error as a result of singularities is greater for larger h and 
this is reflected in these increasing corrections for higher h. 

7. Summary 

The most impressive application of the modified Borel resummation technique in this 
paper is not to look for an approximate value of a resummed expansion but instead 
to use the analytic properties of a function to tune for the correct boundary condition 
(section H]). This method was ffist outlined in p] and in this paper we have given 
numerical evidence to explain why the technique works. Due to the arbitrary levels of 
accuracy that may be obtained through this method it is worthy of futher study. It has 
already been applied to anharmonic oscillators with higher order potentials however we 
believe it may be adaptable to solve other types of problems. 
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The concept of direct resummation is also important. We have used it to 
numerically investigate the location of zeros in solutions to the anharmonic differential 
equation and therefore explain how the tuning technique selects the solution with the 
correct boundary condition. We found that when the correct boundary condition is 
not satisfied the logarithm of the wavefunction contains singularities either on the real 
axis or in the complex plane depending on whether the energy chosen is too large or 
too small. The existence of these singularities causes an easily observable contribution 
though the modified Borel resummation procedure in the form of a rapidly increasing 
or decreasing curve. 

In other physical problems we are not always in a position to apply a tuning 
type method where we look for the correct large x behaviour. We have shown 
that the modified Borel summation technique can produce reasonably accurate results 
from perturbative expansions that are ordinarily divergent. This has been shown 
by considering two different perturbative expansions of the wavefunction. The two 
expansions have different singularity locations and residues but in both cases these 
are the half x plane with negative real part. The semi classical expansion in h was 
particularly useful because it separates the g^^^ type behaviour for large g from the 
expansion we intend to resum. Therefore we get much better results for this part of the 
spectrum than if we used the Bender Wu coupling expansion. 

We have taken the modified Borel technique and improved it by comparing two 
different methods for approximating the contour integral. Both techniques have their 
advantages depending on the location of any singularity contributions and how changing 
the a parameter causes them to move through the complex plane. We have investigated 
the major sources of error and improved the accuracy of the procedure by subtracting 
the leading order corrections. Also, in contrast to standard Borel summation we do not 
require an analytic continuation of the Borel sum since this is encoded in the contour 
integral. 
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